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Abstract 

The field theory of self-avoiding tethered membranes still poses major challenges. In this 
article, we report progress on the toy-model of a manifold repelled by a single point. Our 
approach allows to sum the perturbation expansion in the strength go of the interaction 
exactly in the limit of internal dimension D — > 2, yielding an analytic solution for the 
strong-coupling limit. This analytic solution is the starting point for an expansion in 2 — D, 
which aims at connecting to the well studied case of polymers (D = 1). We give results 
to fourth order in 2 — D, where the dependence on g is again summed exactly. As an 
application, we discuss plaquette density functions, and propose a Monte-Carlo experiment 
to test our results. These methods should also allow to shed light on the more complex 
problem of self-avoiding manifolds. 
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1 Introduction 

One major problem in statistical physics is the effect of interactions on the thermodynamical 
properties of extended fluctuating geometric objects. In general, multi-particle attractive or re- 
pulsive interactions are involved. One may divide these into two classes: Either (i) one may 
study the interaction of a single fluctuating object with itself as for instance the well known 
excluded volume interaction between any two monomers in a long polymer chain in a good sol- 
vent. There, the interaction leads to universal long-distance properties of chains as for example 
the anomalous scaling of the mean squared end-to-end distance. Or (ii), the interaction may act 
between different manifolds or between a single manifold and a fixed non-fluctuating object. It 
is then interesting to study how thermal fluctuations affect the depinning of the manifold from 
an attractive substrate as well as the steric repulsions from a wall. Finally, both cases can appear 
together. 

Whatever the situation is, it is usually well understood as long as the fluctuating objects are one- 
dimensional Referring to the example mentioned above, the long-distance properties of 
self-avoiding polymers can be analyzed with renormalization group techniques [5-7], either in 
the continuous Edwards Hamiltonian [8], 

H[r\ = l - J (Vr») 2 + | J J S d (f(x) - f(y)) , (1.1) 

xeM xeM yeM 

or by mapping this model on a local O(N) symmetric (/? 4 -theory in the limit of N = compo- 
nents [9,3,1]. The critical exponents describing the long-distance properties are related to the 
critical exponents of the corresponding N- vector model at the critical point. What makes (1.1) 
a non-standard theory is that the interaction is non-local, and not a polynomial of the field. 
Obtaining the corresponding results for membranes poses considerable challenges. The general- 
ization of polymers to 2D-surfaces are crystalline fixed-connectivity membranes as they appear 
for instance in the spectrin network of cell membranes. Considering "phantom" membranes 
which can freely fold into itself, the existence of a bending rigidity induced phase transition 
separating a high rigidity, low temperature flat phase from a low rigidity, high temperature 
crumpled phase is well established [10-15]. This is in contrast to polymers, which are always 
crumpled on large scales. The scaling properties of the crumpled phase of phantom membranes 
are described by the 2D generalization of the free field part in (1.1). Taking self-avoidance into 
account, which is modeled in (1.1) through the short-range two-body interaction, we expect 
more swollen manifolds than those predicted by the free theory, which will be expressed in a 
non-trivial radius of gyration exponent v: 

R g ~ L u , < v < 1 , (1.2) 

where L denotes the linear internal size of the membrane, and the radius of gyration R g is 
obtained from the effective extend of the membrane in external space. In the case of polymers 
R g scales like the end-to-end distance. Much effort has been spent on calculating corrections 
to the radius of gyration exponent within an expansion in the deviation e from the critical space 
dimension [16,17]. These calculations can not be performed directly for membrane-dimension 
D = 2, since the naive scaling dimension of the coupling in (1.1) equals 

e(D,d):=[b ]=2D-^^d, (1.3) 
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where d denotes the dimension of the embedding space, such that e is always non-zero as 
D — > 2, for any embedding dimension d. Equivalently, the critical embedding dimension 
defined through e(D, d c (D)) = becomes infinity in this limit. The reason is that the non- 
self-avoiding membrane densely fills out the embedding space, such that it always "sees" the 
interaction. A way to circumvent this problem is to set up the expansion about any point (D < 
2,d c (D)) and to extrapolate along an appropriate path in the d) -plane to the physically 
interesting point (D,d) = (2,3) [18-24]. To second order in e one then finds a radius of 
gyration exponent of v ps 0.86 [16,17] , which is a strong correction with respect to the only 
logarithmic dependence in the non-interacting theory, and indicates the existence of a crumpled 
phase, for which v > 2/3 follows from the fact that a membrane has a finite volume. 
However, there is no evidence for a crumpled phase in experiments [25-28] . Latest Monte-Carlo 
simulations on plaquette-models [29-32] starting from a discretization of the 2D generalized 
Hamiltonian (1.1) with system sizes of up to ps 17000 plaquettes shows considerable evidence 
for a vanishing of the above mentioned crumpling transition in the presence of self-avoidance, 
such that even on large scales fixed-connectivity membranes stay always flat with a radius of 
gyration exponent of v pa 1. The final purpose of this work is to develop techniques which 
allow to go beyond the two-loop result. So far, we developed such techniques for a simplified 
model, which reduces the non-local self-avoiding interaction in (1.1) to self-avoidance with 
only a single point, e.g. the origin o in the membrane: 



This is a special case of a more general model of the interaction of a phantom tethered mani- 
fold with a single point in embedding space, which is related to case (ii). The corresponding 
physical situation to think of is the binding and unbinding of a long chain as e.g. a polymer or a 
membrane from a wall or the wetting of an interface. More precisely, we study the interaction 
of a single freely fluctuating manifold with another non-fluctuating, fixed object. Depending 
on whether the interaction is attractive or repulsive, one can distinguish two different scenarios: 
One may either observe a derealization transition from an attractive substrate as in wetting 
phenomena or steric repulsions by a fluctuating manifold. Both cases have in common that ex- 
cluded volume effects become important. These scenarios have already been discussed in [33]. 
The result of [33] is the complete resummation of the perturbation series for the effective cou- 
pling in the case of 2_D-membranes. The long-distance behavior of the resummed theory turned 
out to be non-trivial in the sense that it emerged from the limiting behavior of a scale invariant 
theory resulting in an effective coupling growing logarithmically instead of approaching some 
finite fixed-point value as one would expect it. This, together with the extremely slow conver- 
gence of the perturbation series makes the analysis of the fully resummed theory a must, because 
all finite loop calculations fail to extract the correct large distance properties. The importance 
of the result becomes evident as soon as one compares it with extrapolations obtained from the 
^-expansion at the 2-loop level [33]. Besides the necessity of calculating diagrams numerically 
with considerably raising effort as the loop order becomes higher, the ^-expansion has turned 
out not to be able to make reliable predictions for D < 2. This problem persisted, though we 
exploited the freedom to set up the expansion about any point (D < 2, d c (D)), d c (D) = 
being the critical embedding dimension for given internal dimension D and to expand both in D 




(1.4) 



xeM 



xeM 
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and d along any appropriate extrapolation path to some physically interesting point (D = 2,d). 
As soon as D approaches 2 the result was always strongly dependent on the selected expansion 
point. 

The aim of this paper is two-fold: First, we reconsider the techniques to perform loop cal- 
culations within a "massive scheme", that is on a manifold of finite size and with fixed space 
dimension < D < 2 and d. We show that the perturbation series for the effective coupling 
can be completely summed in D = 2 and analyze the long-distance properties in this limit. In 
addition to [33], instead of analytically continuing loop integrals to D = 2 from below we also 
perform calculations directly in 2D, which need an explicit short-distance (UV)-cutoff. It turns 
out that results in D = 2 are independent of the procedure, i.e. they are universal. 
Second, we construct a systematic expansion of the effective coupling in powers of 2 — D. 
Such an expansion is based on our techniques to resum the perturbation series at each order in 
2 — D. A first attempt to go beyond D = 2 has already been made in [33]. However, there, the 
effect of the boundaries of the finite manifolds was not taken properly into account, a problem 
that has now been circumvented by considering closed manifolds. We specialize to a toroidal 
internal topology corresponding to periodic boundary conditions. Of course, the propagator of 
the perturbation series needs to be modified, and diagrams become more difficult to calculate. 
Slightly below D = 2 we expect power-law behavior of the effective coupling. We present a 
possible ansatz for the exact effective coupling as a function of the internal dimension D < 2, 
which is consistent with the expansion m2 — D. However, it remains an open problem to obtain 
more information about the power-law behavior, to make this expansion unique. 
A short account of this work has already appeared in [34]. 

2 Model and physical observables 
2.1 The model 

The problem of a membrane avoiding only a single point (1.4) may at first sight appear artificial. 
Besides of being a toy-model for the analysis of the more interesting case of full self-avoidance 
as discussed above, let us point out that it is, too, a special case of another more general problem, 
which is interesting by its own. Consider a phantom tethered membrane interacting with some 
5-potential located at the origin of the configuration space: The Hamiltonian is given by (0 < 
D < 2) 

W\ = \ f (Vr») 2 + <?o J S d (f(x)), (2.1) 

where any point in the membrane is labeled by some D-component vector x, and its position in 
external space is given by the rf-component field f(x), 

r:xeR D — ► f(x) G R d . (2.2) 

The partition function is defined as 




(2.3) 
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To remove the translational 0-mode, we will consider 



= jv[r\ 5(r(y)) exp(-W[f|). 



(2.4) 



We now discuss (2.1) in more detail: The first term is the elastic energy of the manifold which 
is entropic in origin. Elasticity and temperature have been scaled to unity. The second term 
models the interaction of the manifold with a single point at the origin in the <i-dimensional 
configurational space. Let us remind [35,33] that the physical interpretation depends on the di- 
mensionality: In the case that IR d is identical to the embedding space, (2.1) describes a phantom 
crumpled manifold interacting with a single defect as sketched in figure 1. However, setting 
d = 1 (2.1) may as well describe a solid-on-solid like fluctuating interface parameterized by 
some displacement field and interacting with a parallel plane (D = 2) as shown in figure 1 . 





Figure 1: Left (a): A D-dimensional manifold (D = 2) interacting with a point in the origin of the configurational 
space M. d . Right (b): A "directed" membrane (interface) interacting with a parallel subspace of same dimension 
D. 

The coupling constant g may either be positive (repulsive interaction) or negative (attractive 
interaction). We now give the dimensional analysis. In internal space units, the engineering 
dimensions are 

dim[x] = 1 

j- r-a 2-D 

v := dim[r] = 



e :- 



dim 



d v x S d (f(x)) 



M 



D-vd. 



The interaction is naively relevant for e > 0, i.e. d < d c with (see figure 2) 

2D 



d n 



2-D 



(2.5) 



(2.6) 



irrelevant for e < and marginal for e — 0. It has been shown [35,36] that the model is 
renormalizable for < D < 2 and e > 0. Results for negative e are obtained via analytical 
continuation. One can define the renormalized coupling g as 
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D 




Figure 2: Critical line defined through e = d c (D) = t^tj ■ The interaction is relevant for points that lie above 
that line. 



9 ■-- 



V M 



[Z(0) - Z(g )] U , 



(2.7) 



where Vm denotes the internal volume of the manifold. The normalization J\f depends on the 
definition of the path-integral (but not on L) and is chosen such that 



g = g L s + O(g 2 ). 

Universal quantities emerge at fixed-points of the /5-function, which is defined as 

dg 



(3(g) := —L 



dL 



(2.8) 



(2.9) 



.<?() 



The /^-function thus describes, how the effective coupling g changes under scale transforma- 
tions, while keeping the bare coupling g fixed. Let us give the 1-loop result, see e.g. [35,36,33]. 



$(g) 



$(g) 

















y. 







g 





g 



(a) 



(b) 



(c) 



Figure 3: RG-function and flow for increasing manifold size L for the dimensionless renormalized coupling g: (a) 
in the case e>0, (b) in the case e<0, (c) in the case e = 0. 
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It reads 



(3(g) = -eg + ±g 2 + 0(g 3 ) 



(2.10) 



where g is the dimensionless renormalized coupling. Apart from the trivial solution, g — 0, the 
flow equation given by (2.9) and (2. 10) has a non-trivial fixed point at the zero of the /3-function 



We shall show below that the scaling behavior is described by the slope of the RG-function 
at the fixed point, which is universal as a consequence of renormalizabilty. The long-distance 
behavior is then governed by the 5-interaction as considered in our model (2.1), which is the 
most relevant operator at large scales. Let us now discuss possible physical situations (see fig. 



(a) e>0: The RG-flow has an infrared stable fixed point at g* > and an IR-unstable fixed 
point at g — 0. The latter describes an unbinding transition whose critical properties 
are given by the non-interacting system, while the non-trivial IR stable fixed point deter- 
mines the long-distance properties of the delocalized state, the long-range repulsive force 
exerted by the fluctuating manifold on the origin - which we remind may be a point, a 
line or a plane. 

(b) s<0: Now, the long-distance behavior is Gaussian, while the unbinding transition occurs 
at some finite value of the attractive potential, g*<0, which corresponds to an infrared 
unstable fixed point of the /5-function. Below g* the RG-flow is to strong coupling and 
the manifold is always attracted. 

(c) e = 0: This is the marginal situation, where the transition takes place at g* = 0; we 
expect logarithmic corrections to scaling. 

These scenarios and possible observables have already been discussed in [33]. Here we want 
to specialize to the case of a membrane avoiding a single point. It turns out that this situation 
allows to calculate observables staying non-singular even for 2D membranes and which are 
accessible to a Monte-Carlo experiment. 



Interesting physical observables for a membrane avoiding a single point are the plaquettes- 
density functions at the repelling point. Generally, these are defined as follows: 



9* = 2e + 0(e 2 ). 



(2.11) 



3): 



2.2 Plaquettes-density correlation functions 




(2.12) 



where the expectation value (-) is taken within the pinned ensemble as defined in (2.4). The 
quantity, which is accessible in perturbation theory, is the effective coupling as defined in (2.7). 
It can be considered as a generating function for observables like (2.12). Let us first show how to 
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obtain the constrained partition function (2.4) from (2.7): Since we consider closed manifolds, 
internal translational invariance implies 



1 d 



<7( \ 1 d ^ L ~ 



(2.13) 



where g is the renormalized or effective coupling defined in (2.7) and Vm denotes the internal 
volume of the membrane. Introducing the dimensionless bare coupling, 



z := g L £ , 

(2.13) can be written in terms of dimensionless quantities as 



(2.14) 



(2.15) 



where J\f has been set to unity. In the same way, all observables of the type (2. 12) can be easily 
derived from g according to: 



L le e fdg 



5a dz 1 V dz 

dz 



(2.16) 



An observable, which is accessible through Monte-Carlo simulations, should be expressable by 
universal quantities like the slope of the renormalization /3-function (2.4). All these quantities 
can be obtained from appropriate derivatives of the effective coupling g with respect to z: The 
/3-function can be written in terms of the bare coupling as 



(Note that in a slight abuse of notation, we write (3{z) - 
fixed point, which is defined as 

d(5{g) 



is obtained from 



u(z) 



dg 

—ez df3(z) 



(2.17) 

f3(g(z)).) The universal slope at the 

(2.18) 



P(z) Oz (Z19) 
in the limit z — > oo. We furthermore need the second derivative of the RG-flow function with 
respect to the effective coupling, which is defined as 

d 2 (3(g) 2 _>oo -ez duo(z) 



dg 2 



(3{z) dz 



(2.20) 



Let us now show that the universal slope (2.18) is accessible through the measurement of certain 
combinations of observables of type (2.12). For this purpose we need the plaquettes-density 
(£ — 1) and the density-density function (£ — 2), which are obtained after some straight forward, 
but tedious algebra from the above definitions: 



« = — ( 1 + 

go 



g \ e) 



9l 



9l 



3u u 
2 + — + - 

e e- 



(2.21) 
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These quantities depend on the bare coupling g , which is not accessible. Instead, consider the 
following ratio: 

(n) z _>o 



o 




(2.22) 



which obviously is universal. 



2.3 Delocalization transition 

For completeness let us shortly discuss the physical situation at the UV-stable fixed point in 
figure 3. The fixed point corresponds to a delocalization transition of the manifold, which is at 
vanishing coupling g*=0 for s>0 and at some finite attractive coupling g*<0 for e<0. 
In the localized phase g<g*, correlation functions such as ([f(x) — r(y)] 2 ) and the associated 
correlation length £y (in the D-dimensional internal space) should be finite, as well as the radius 
of gyration Approaching the transition these quantities diverge as [37] 

£H ~ (g*-g)-^ , a ~ (9*-g)- u± ■ (2.23) 

Since £j_ ~ the exponents v\\ and v±_ are related through 

v± = , (2.24) 

v being the dimension of the field (2.5). 

Furthermore, they are related to the correction-to-scaling exponent oo: 

v \\ = , "± = r^T • (2-25) 

Note that u(g*)<0 at the transition. Specializing to (D, d) = (l, 1), we find 

u ± = 1 , = 2. (2.26) 

These exponents are also valid for the delocalization transition of a 1 -dimensional interface 
from an attractive hard wall in 2-dimensional bulk space [38,37,39,33]. 



3 Complete summation of the perturbation series 

3.1 Perturbation theory 

In (2.2) we saw that physical observables can be derived from the renormalized coupling g (2.7). 
To obtain g we need the perturbation series of the partition function Z (2.3): 

where ^ 

Zn = ( II / * d (Kzi))) , ^>0, (3.2) 
\ i=i Jxi /o 
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and the normalization of the 5-distribution has been chosen to be 

5 d {r(x)) := (Au) d l 2 5{r(x)) = I e lkr ^ (3.3) 

J h 

with 

J := n'^ 2 J d d k . (3.4) 
The advantage of these normalizations is that 



[e- k2 = l. 
Jk 



(3.5) 



ik 

Accordingly, the perturbation expansion of the effective coupling (2.7) reads 

\N l N + 1 





Performing the averages within the Gaussian theory with normalization 

1 



Vm J x \ I o 

one arrives at 

'N+l „ „\ N+l 



5 d (r(x))) =1, (3.7) 



Vm ^ (N+l)\ 

where 



N=0 v > \ i=l t " / i=l 



C(x 4 -a;,):=l<(r(^)-r(^)) 2 > (3.9) 

denotes the correlator, and the 5 d (J2i stems from the integration over the global translation. 
Shifting 

N 

^n+i k N+ i — hi , (3.10) 



=i 



the quadratic form in (3.8) transforms to 

N+l 

2 



- ^ kikjC(xi Xj) 

j=l i,j=l 

(3.11) 

Integrating out the momenta k 1: . . . , A:at + i in (3.8), one obtains 

^ )=, £SS(n/)< drt » r * > - (3 - l2) 
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where we have factored out U from the loop integration (such that the integrals now run over a 
torus of size 1), and the matrix elements 2)^ are 

®ij = ^[C{x N+1 -x i ) J r C(x N+1 -x j )-C(x i -x j )\ . (3.13) 



3.2 Complete summation in fixed internal space dimension D = 2 

Let us compute the iV-loop order of (3.12): The behavior of the propagator C(x) for arguments 
x large compared to a is of the form 

C(x)=c + ^\n-, (3.14) 

where c denotes some positive constant (note C (x) > 0), and the logarithmic growth (for large 
x) is universal (see appendix A). In D = 2 we need an additional short-distance cutoff a, which 
we want to take to 0. We can (somehow arbitrary) decompose 

N 

det£> = (JJS^detS) . (3.15) 

i=i 

In the limit of a -> each C(x) = j- \n(L/a) + 0(a°), such that 



1 | C(x N+1 -Xj) - C(xi-Xj) 



C^jv+i-^i) 

= 1 ■ (3.16) 

(det D)" d/2 =: 7jv(L/a) = /f (L/a) (det &Vy d/2 ■ (3.17) 

The matrix 2) (0) denotes the limit a — > of (3.16). It can be written as I) (0) = |(I + ATP), 
where I denotes the identity and P the projector onto (1,1,. ..,1), whose image has dimension 

1, such that det = ^ff- [33]. Furthermore, to one loop I^L/a) a =° Ci(ln ^)~ d / 2 , where c x 
denotes some (finite) constant. One then arrives at 




A7i(I^v)d/2+i • 



g( z ) = ^ v u^i • (3 - 18) 



iV=0 

A factor Ci2 d / 2 has been absorbed into a rescaling of both z and g. 



3.3 Asymptotic scaling behavior 

In the following we will analyze the limit of large z (strong repulsion), which also is the scaling 
behavior of infinitely large membranes. We need an analytical expression for sums like (3.18) 
in the limit of large z. Later, it will turn out that allowing for small deviations 2 — D > only 
slightly more general sums will arise. 
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We claim that for all k, d > 



E 

N=0 



N\(k + N) d / 2 r(±)Jo 



1 f°° 

/ dr r d/2 ~ l exp[-z e~ r - kr] . 
I Jo 



(3.19) 



This can be proven as follows 
1 



r(f) 



poo i / \Af poo 

/ dr r d ' 2 - 1 exp(-z e~ r - kr) = V / dr r «*/2-i e -("+*)'- 

CO 



This integral-representation is not the most practical for our purpose. It is better to set r — > s : = 
e~ r which yields 



E 

N=0 



N\(k + N) d / 2 



r(|) 7o 1 ; 



(3.20) 



This formula is already very useful for some purposes. It is still advantageous to make a second 
variable-transformation s — > y := sz, yielding 



E 

AT=0 



N\(k + N) d / 2 



(lnz)^ 2 - 1 



k-1 



1 _ ln ^ ^ 

lnz/ 



d/2-l 



Finally we remark that we usually have the following combination 



(lnz 



,d/2-l 



AM(A; + iV) d / 2 



fc-i 



1 - 



lnj/ 



d/2-l 



(3.21) 



(3.22) 



JV=0 - ■ - v~ ' - ■ / v 2 7 

It satisfies the following simple recursion relation, which is helpful to calculate the /5-function: 



(3.23) 



The derivative above can also be rewritten as 

d 



(3.24) 

such that one obtains a useful formula in order to isolate the dominant behavior for large z: 

ft +1 (z) = kf d (z) - ft 2 (z) . (3.25) 
From (3.19) f^(z) > for all k, d > and the behavior for large z is obtained by expanding 
^ small £ 



(lnz 



id/2-1 



dyy 



1 /d 



lnz V2 



- 1 



+ 



1 



(lnz)' 



dy y k 1 In y e v 



(3.26) 
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The result is 



m = 



(In z) d/2 ~ 1 T(k) 



r(|) 



1 d-2T'(k) 



+ ... 



In; 2 T(k) 

With the above notations, the sum (3.18) expressing g as a function of z becomes 



(3.27) 



g(z,L/a) 



T \ d/2 

In-) ft +2 
a J 



In 



-d/2 



(3.28) 



in the limit D = 2. 

It is now easy to analyze the long-distance behavior in this limit. First, we observe that accord- 
ing to (3.27) the effective coupling diverges logarithmically for all external dimensions d > 0: 



g(z,L/a) 



(In f) 



d/2 



ln(,(ln^ /2 ); 



d/2 



(3.29) 



This is in contrast to the one-loop result as stated in (2. 11), which is exact for polymers (D = 1) 
and which stays qualitatively valid as long as D < 2. This follows from the renormalizability 
of the theory [35] for sufficiently small e > 0. A finite limit g(z — > oo) = g* signals a scale 
invariant theory. In (3.29) we have found the limiting behavior of the latter. Consequently, we 
expect the correction-to-scaling exponent uj to be always zero in D = 2. In order to check that 
let us first compute the renormalization (3-f unction in terms of the bare coupling as in (2.17), 
which can be immediately derived with the help of relation (3.23) 1 : 



1 



r(|) 



[Hz')] 



d/2-l 



(3.30) 



where we have introduced rescaled couplings z' := g L £ '(In ^) d l 2 and g' — g (In ^) d l 2 . Its 
derivative with respect to the renormalized coupling is found as a function of the bare coupling 
(2.19) to be 



uj(z') 



-ez' d(3{z') 
/3(z') dz' 



—e- 



y< -A 

" dz 



2-d 
2 In(z') 



. 



(3.31) 



Note that the qualitative behavior of the /5-function changes depending on the external dimen- 
sion d, approaching asymptotically zero below d = 2 and being divergent above. 
In the limit of large bare couplings one may as well give the RG-function in terms of the ef- 
fective (renormalized) coupling simply by inverting the asymptotic expression in (3.29) and 
inserting it into (3.30), with the result: 



— >oo 

~ — e 



(r(d±2))i-2/ d 
r(|) 



l-2/<2 



(3.32) 



'Note that our definition (3(z') = —ez'dg'/dz' is strictly speaking equivalent to definining the /3-function as 
f3(g) := (— — fflgjj) | g, instead of (2.9). (Note that the derivative w.r.t. a disappears for D < 2.) The 



natural combination in D = 2 is z' = g^L e (in ^) instead of z = goL £ , and normalizations such that 
g'(z') = z' + 0(z' 2 ) does not explicitly depend on L or a. The chosen defintions avoid unnecessary technical 
complications, but do not change the physics of the problem. 
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Figure 4: (3— function in terms of the renormalized coupling g truncated at order 160, Pade-resummed, and plotting 
only that part for which the truncated series converges. (This can e.g. be tested by taking away the last few terms of 
the series.) This is compared to the asymptotic behavior (3.32) (proportional to 1/g for large g). d is set to 1, and 
we used the diagonal (80,80)-Pade approximant, which was find to converge best. (The non-resummed expression 
starts to diverge already at g as 1.8 at this order.) 

It is interesting to compare the true asymptotic behavior of the completely resummed pertur- 
bation series as found above with predictions taking only finite loop orders into account: If 
one tries to invert (3.18) and truncates it at some finite order, it is at least possible to reach 
the asymptotic regime (3.32) - however, for large g the truncated (5— function does not con- 
verge to the true /3-function and thus strongly deviates from the true behavior. In figure 4 the 
Pade-resummed truncated /3-function up to order g 160 in d = 1 is compared with the asymp- 
totic flow-function. One notices that the truncated /3-function even though improved through a 
Pade-Resummation hardly gets into touch with the asymptotic regime. The same applies to the 
slope-function co(g), which is not shown in figure 4. Let us finally state the expected behavior 
of the plaquettes-density functions in the limit of large membranes. For the plaquettes-density 
at the repelling fixed-point we find in this limit: 



Note that in the absence of the repelling interaction this quantity would diverge in this limit. 
This follows from dimensional grounds, since then 



In (3.33) we found the largest possible depopulation of monomers at the defect potential in the 
case of a relevant interaction (e > 0). As we discussed in (2.2) a measurable quantity should be 
the following ratio (2.22), which in the case of 2D-membranes becomes in the limit z — > oo: 




(3.33) 



(n) ~ U . 



(3.34) 




(3.35) 
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Figure 5: Qualitative behavior for the /3-function in D = 1, D = 2 and result anticipated for D w 1.5. 



which can be compared with the 1-loop prediction (which is exact for polymers): 



rr 



(1-loop) 



(3.36) 



4 Crossover to polymers 

Let us now analyze the theory below D = 2. Due to the renormalizability in < D < 2 and the 
existence of an ^-expansion we expect the renormalized coupling to reach a finite fixed point 
in the strong coupling limit as soon as D < 2. This approach is characterized by a power-law 
decay of the form 

g (z) =g* + S(\nz) z~"' e + 0{z~^ £ ) , (4.1) 

where S is some scaling-function growing at most sub-exponentially and uoi > uo > 0, with uo 
defined in (2.18). 

Our ultimate aim is to draw information from an expansion in powers of 2 — D of the effective 
coupling about the correction-to-scaling exponent uj in (4.1) for D < 2. The scale invariant 
behavior below D = 2 results in a finite fixed point of the renormalization /3-function as a func- 
tion of the effective coupling. The qualitative behavior of the /5-function is sketched in fig. (5). 



4.1 (2-£>)-expansion on the torus 

In order to gain information about g below D = 2 one has to expand the loop integrand 
(det D)~ d / 2 (3.12) in powers of 2— D. For convenience, we take a — > 0. The propagator takes 
in infinite D-space the form C(x) = \x\ 2 ~ D / {S D {2-D)), where S D = 2vr D / 2 /r(f ) denotes 
the volume of the D-dimensional unit-sphere. The factor (S£>(2— D))^ 1 replaces ln(^) and is 
absorbed into a rescaling of the field and the coupling according to r — > r (S D (2—D)) 1 / 2 and 
9o -> (^d (2-D)) d / 2 , such that the factors of (In ^)- d / 2 in (3.18) and (3.29) disappear. The 
propagator in the rescaled variable can then be written as 

C{x) = 1 + (2 - D) C{x) , (4.2) 

where for convenience of notation we allow C(x) to depend itself on D. 
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Of course, on a closed manifold of finite size, C(x) is modified, but the form (4.2) is inde- 
pendent of the shape of the manifold. Accordingly, one may expand the matrix D as 



D = £ {0) + (2-D) D , 



(4.3) 



where is defined as before and coincides with the limit D^2 when inserting the above 
C(x) into D. Moreover, D is of the same form as D, but each C(x) has been replaced with 

C(ar): 



Then, 



% = -[C(x N+ i - Xi) + C(x N+1 - Xj) + C(xi - Xj)\ . 



detD = det£ (0) exp jlY ln(l + (2 - D)^^)]" 1 ©)] } , 



where [D (0) ] _1 = 2(I-^-P) denotes the inverse matrix of !D (0) . 
Denoting 

M := [D^]- 1 ® , 
we expand the determinant in (4.5) up to fourth order in 2—D: 



(4.4) 
(4.5) 

(4.6) 



1 - ^ [(2 - D)TvWl 
TrfH 2 + v "~" ; Tr9rt 3 - 



[det(D)]~ d/2 = [det{T) i0) )]- d/2 
(2 , (2-^ 3 



(2-DY 



TrOrt 4 



rf 2 
48 



(2 - £)) Tr Ort - (2 - -D) Tr9rt Tr9rt + (2 - D) 1 



-Tr 2 9rt 2 + -Tr9rtTr9rt 3 

4 3 



(2 - D) 3 Tr 3 9rt - (2 - £>) 4 -Tr 2 Tt TrOrt 2 



+ 3^(2-£) 4 Tr 4 9rt 



0((2-D) 5 ). 

(4.7) 



The first step in the analysis will be to obtain the resummed perturbation series of the effective 
coupling up to fourth order in 2 — D. That is, we have to insert (4.7) into (4.5), calculate the 
corresponding loop integrals at each order of perturbation theory, insert the result into (3.12) 
and sum the appearing series to all orders. 

Let us start with the first-order term in 2 - D from (4.7). We only need 9rt = [D^ 1 !), 
which reads 



N 



= (p)' ']- 1 ©).. = ( 22),, - — $> fc ) L 



-2u 



k=l 



The trace of (4.8) can easily be performed, with the result 

n 2 N N 

T 



TrOrt 



i=l i=l fc=l 



2i/ 



(4.8) 



(4.9) 
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Figure 6: Regularization scheme for the iV-loop diagrams on manifolds with toroidal topology (periodic boundary 
conditions). Here: D = 2 

In each order of perturbation theory we have to integrate the expression (4.7) over internal 
distances. These integrals have to be regularized in the infrared through an appropriate IR cut- 
off. We are considering a finite manifold of toroidal topology. The precise form of the correlator 
on the torus will only later enter into the calculation. 

To simplify the calculations, we further introduce the following notation: 



f(x il ,...,x ik ) := / •••/ f(x h ,...,x ik ) (4.10) 
with the internal integrations defined as 

/= L £ / , / := integral over the torus with L = 1 , (4.11) 
J X J X 

such that the overbar in (4.10) can be thought of as an averaging procedure, and especially 

1 = 1. (4.12) 

Thanks to our regularization prescription the integral of (4.9) over internal points can be re- 
placed by L ND (for the integration measure) times 

2N 2 



TrSTJl = — — C(x N+1 - Xi) 



2N N(N-l) 

= tTn c{xn+1 ~ Xi) + i + N c{Xl ~ Xj) ■ (4 ' 13) 
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Due to the internal symmetry of the closed manifolds which we consider the expression above 
can be further simplified, since 



C(x N+1 - Xi 
Introducing a diagrammatic notation 



[ Xi Xj 



C(x) . 



(4.14) 




C(x), 



(4.15) 



the iV-loop integral reads up to first order m2—D 



I 



(det D) 



-d/2 = ^-Ne 



1 + N 
2 N 



-d/2 
X 



l--(2-D)[N v ^ 




+ 0{{2-D) 2 ) 



. (4.16) 



For the further analysis we will not only need (4.13), but also the terms appearing to higher 
order in 2 — D in (4.7). We derived expressions like (4.13) for Tr 2 M and TrM 2 and all terms 
up to fourth order in 2 — D with a Mathematica®-program. It is based on the fact that all terms 
to appear in the expansion (4.7) are of the form Tr n M m or products of the latter and therefore 
can be written as P(N)/(N+l) k , where n,m,k G N and P(iV) is some polynomial in N. It 
will turn out soon that it is convenient to expand the polynomial P(iV) in terms of the following 
base: 



We obtain: 



{1, N, iV(iV-l), N(N-l)(N-2), . . . , H(N-j), ...} 

3=0 



-2N(N-1) - N(N-l)(N-2) 



(4.17) 



TrOJl 2 = 



1+N 

2N + 3N(N-1) + N(N-l)(N-2) 
T+N 



C 2 (x) 



and 



Tr 2 9Jl 



4N(N-1) + N(N-l)(N-2) 
T+N 



C(x) + 



2N 
1+N 



C 2 (x) 



(4.18) 



(4.19) 



Diagrammatically, the averages can be rewritten as 




and 



:=C(x) , 



(4.20) 



(4.21) 



Like in the case of the first order diagram (4.18) and (4.19) are highly simplified as compared 
to an open manifold, see our treatment in [33]. 
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Let us shortly discuss the reason for (4.17): Inserting (4.18) and (4.19) into the perturbation 
series and summing all loop orders, the following series types will appear: 



I]! !j(N-i)(-z) N = k k ^ (-z) N = k d+2{j+1) 

\(N+iyi 2+ i +i n N K N + k ) d/2+j+1 ~ 



Z ^ .Y! 

N=0 



(4.22) 



We may therefore identify the resummed series with a function that we know already fairly 
well, in particular we know its strong coupling behavior. It is furthermore convenient to reduce 
all functions ft>\^ +1 \ z ) to sums of functions ff +2 ^ +1 \z) exploiting the formula (3.25). 

4.2 Resummed contributions to the expansion in 2 - D up to fourth order 

We are now almost in the position to state all resummed contributions up to fourth order in 
2 — D. Let us first state all necessary diagrams: 




which contributes to first order in 2 — D. To second order one needs in addition 




To third order diagrams with new topology are 




C^(xj Xj) , 




= Cixi-x^Cixj-x^Cixi-Xk) 



Finally, to fourth order arise: 





^(xi-x^Cixk-x^Cixk-Xi) , 



C(x i -x i )C(x fc -a; i )C(x i -x fc )C(a; i -a; i ) . 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



If one calculates diagrams, it will turn out that it is to some extend more convenient to express 
the above averages in terms of averages over a connected correlation function, which is defined 
as 

C c (x) := C(x) - C , (4.27) 

such that for instance 

Cf = - C 2 . (4.28) 
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Furthermore, we will need: 



c 

and 



C3 = C 3 -3CC 2 + 2C 3 (4.29) 



= C(xj— Xj)C(xj— Xfe)C(xjfc— Xj) + 3C C(xi-Xj) — 3CC(xi—Xj)C(xj—Xk) — C 
= - C 3 , (4.30) 

where Xi, Xj , Xk are distinct points, and the average is over their positions. In (4.30) we exploited 
the symmetry of the closed manifold, and the definition of C A is self-evident. Furthermore, we 
will need to fourth order in 2 — D: 

C* = C* + 12C2C 2 -4C3C-3C 4 , (4.31) 



Q = CciXi-X^CciXj-X^CciXk-X^CdXi-Xi) 

= + 5C 4 (4.32) 

and 



C/f- = C2(x i -x i )C c (xi-a: A: )C c (a; fe -a; i ) 



= C°> - 2C A C - C 2 C + 2C . (4.33) 

Let us now state all terms which appear in the expansion of the renormalized coupling g(z) up 
to fourth order in 2 — D according to (4.7). We have to calculate at order N of perturbation 
theory: 



Tr = iVC. (4.34) 

Inserting this into the perturbation series and summing up the resulting terms to all orders in iV 
generates the following contributions in the 2 — D-expansion of the renormalized coupling: 



E(detgW)^ ^?^ =Cf*»(z) -Cf*(z) , (4.35) 



(N+l) 

N=l v 1 

which contributes to first order in 2 — D. 

To second order in 2 — D, we have (4.18) providing 



^(detg'°')-/ 2 lt -°" ■ J> = 2Cf + (-4Cf + T*)f?\z) 



and (4.19) providing 

\N+1 



+( _C 2 + 3C 2 ) f*{z) - C 2 ft 2 (z) , (4.36) 



^(detD(°))'^ Tr + = - (2q +cvr(^) 

+2C 2 / 1 d (z)-C 2 /i d - 2 (^) • (4-37) 
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Let us now state the terms at third order in 2 — D, which we derived with the help of a 
Mathematica®- program (N is the loop order.): 



+(-10C| + 36C C A + 6C C*)ft +2 (z) 

+(9Cf - 32C? - 12CCf + C 3 )f?(z) 

+(3Cf - 17C? - 9C Cf + C 3 )/i d_2 (^) 

-3(2C? + C Cf)A d - 4 (z) + C? ft\z) , (4.38) 



AT=1 



(iV+1)! 

+ (-8C3 + 32Cf + 2CC2)/f +4 (z) 
+ (6Cf - 20C? + 2C Cf - 6C 3 ) f? +2 (z) 
+ (-2Cf + 4C? - 7CCf + 2C 3 )/ 1 d (z) 
-(C 3 - 4CCf)/ 1 d - 2 (z) - C Cf /f -4 ^) , (4.39) 



f>t aW)- ^ 3 - 4(Cf - 4C c A )/r (,) 



JV=1 



+(-4C3 + 24C A - 6CC2)/f 



+(-8C A + 12CC2 + C')f« +2 {z) 
-3(C 3 + 2CCf)/f(2) 
+3C 3 / 1 d - 2 (^-C 3 / 1 rf - 4 (^. 



(4.40) 



To fourth order in 2 — D we obtain: 



N=l 



„^r2 



8(222C' + 6C'C2 + 3Cf - C 4 C + 24C C ° > - 36C*)/ 1 d+8 (z) 



,^4 



+4(804C' + 12CC 2 . + 3Cf - 2C^ + 72C C °> - 132Q -ACC 3 C + 16CC* )f^(z 



fci+6/ 



^4 



,^2 



-4(432C - 3CC 2 . - 6Cf + 24C c a> - 72Q - 8CC;? + 40CCf )ft + \z) 
+ (287C 4 



36C C2 



12Cf - 48Q - 16CC3 + 128CCf )/ 1 d+2 (^) 



+ (4C" + 36C*q - 32CCf 
+ (-6C 4 - 2C 2 Cf)/f- 2 (^) + ^ft\z) - ~^ft\z) 



(4.41) 
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S^n ^nm d/2 Tr9Jl 2 Tr 2 <m(-z) N+1 



JV=1 



8(222^ + 6C 2 Cf + 3Cf 2 - Cf + 24C^ - 36Q)/f 
+4(960^ + 24C 2 Cf + cf - 4Cf + 100C^ - 156Q)A d+6 (z) 
-4(714C 4 + 20C 2 Cf + 8Cf 2 - 3Cf + 70Cf r - 116Cf - 4CCf + 12CC?)/ 1 d+4 (z) 
+(889C 4 + 36C 2 Cf - 2Cf 2 - 28C Cf - 4Cf + 80C^ - 144Q + 88C C?)f? +2 (z) 
+(-99C 4 + 3C 2 Cf + 8Cf 2 + 16C Cf - 8C<F + 16Cf - 48C C?)f?(z) 
+(3C 4 - HC 2 Cf - 2Cf 2 - 4C Cf + 8C £^)ff~ 2 {z) 

+(-C 4 - 5C 2 C|)/^ 4 (^) - C 2 C? c ft 6 (z) (4.42) 



•d+8, 



= 8(222^ + 6C 2 Cf + 3Cf 2 - Cf + 24C^ - 36Q)/f+ 8 (z) 

+4(1116C 4 + 36C 2 Cf + 23Cf 2 - 6C| + 128C^" - 180Q + 4C Cf - 16CC?),/f +6 (,z) 
-4(1080C 4 + 47C 2 Cf + 28Cf 2 - 8C| + 132^ - 172Q + 8C Cf - 32C €^)ff +i {z) 
+(2111C 4 + 148C 2 Cf + 44Cf 2 + 24CCf - 24Cf + 272C^ - 328Q - %QC£^)ff +2 {z) 
+(-538C 4 - 74C 2 Cf + lOCf 2 - 8C Cf + IOC* - 72^ + 80Cf + 8C C?)f?(z) 
+(59C 4 + 20C 2 Cf - 15Cf 2 - 2C| + 8C°^ - 8Cf)/^ 2 (^) 

+(-2C 2 Cf + 6Cf )/i d-4 (*) - C?ft 6 (z) (4-43) 



= 8(222^ + 6C 2 Cf + 3Cf - Cf + 24C^ - 36Q) ff +8 (z) 

+4(1038C 4 + 30C 2 Cf + 18Cf 2 - 5C| + 114C^ - 168Q + 2C Cf - 8C C?)f? +6 (z) 
-2(1818^ + 54C 2 Cf + 39Cf 2 + C Cf - 9C| + 204C^ - 294Q - 22C £^)ff +A {z) 



+(1583C 4 + 48C 2 C2 + 36Cf - - 6C* + 186C°> - 258Q + 8CCf )ff +2 {z) 



+(-358C^ - 21C C 2 C - 6C2 + 6CCj? - 48^ + 60Cf - 37CCf )ff(z) 
+(35^ + 12C 2 Cf - 3C Cf + 6C^ - 6Cf + 23C C? ) /f" 2 ^) 

+(-3C 2 - 7CC?)/ 1 d " 4 (2) + CcF/f -6 ^) (4.44) 
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= 8(222C 4 + 6C 2 Cf + 3Cf 2 - Cf + 24C^ - 36C|) f? +8 {z) 

+4(1116C 4 + 36C 2 Cf + 23Cf 2 - 6C* + 128CF - 180Q + 4C Cf - 16CC?)/f +6 (2) 
-4(1110C 4 + 39C 2 Cf + 31Cf 2 + IOC Cf - 7C| + 136C^" - 178Cf - 36CC?)/f +4 (^) 
+(2473^ + 72C 2 Cf + 82Cf + 36C Cf - 16C| + 304C°^ - 396C| - 128C €^)ff +2 {z) 
+(-955^ - 12C 2 Cf - 36Cf 2 - 12C Cf + 5Cf - 92C°^ + 154Cf + 68C C?)f?(z) 
+(288C 4 + 12Cf - Cf + 12C^ - 47C| - 24C Ct)f^ 2 {z) 

+(-60C 4 - 2Cf + lOCf + 4C £^)ff-\z) + (6C 4 - £»)ff-\z) . (4.45) 

4.3 Renormalized coupling 

Combining (3.12), (4.7) and the results (4.35)-(4.45) from the preceeding subsection we may 
now give the exact renormalized coupling to fourth order in 2 — D. For the sake of compactness, 
we introduce a new notation: Since all series contributions are of the form as stated in (4.22), 
we introduce vectors M such that 



n=i UV+ij! j=min y ni n 2 ... ni ) 

= M j , , .f^iz), (4.46) 

/mi m.2 ■ ■ ■ mA max \ </ 1 V / * V / 

\ni n2 ■ ■ ■ ni \ min y 

where max and min are some integers, and summation over the index j is implicit. Inserting the 
results for the resummed series contributions into (4.7) we find for the renormalized coupling 
to fourth order in 2 — D: 

g(z) = ft + \z) - (2 - D) ±\e^ft+*(z) 



+( 2 - of {^w 2 yr j (z) + j m| 2 yt +2j (z)) 



-(2 - Df ( ±m 2a) ft +23 (z) + |m| a^(z) + |V 



The vector entries MP are to be be taken from subsection 4.2. 

It is more convenient to discuss instead of g(z) an integral transform. From the expansion 
of /*(*). namely 



poo 

— drr d/2+j - 1 exp[-ze-' r -r}, (4.48) 
2) Jo 
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and the structure of the expansion of g(z) in powers of 2 — D and the integral representation of 
the ff +2j it follows that the exact renormalized coupling can be written as 



POO 

g(z) = g(D, z) = z dr g(r) exp[- 
Jo 

where g{r) is of the form 



■ZQ 



(4.49) 



g(r) 



„d/2 



1 oo n 



n=0 j=-n n 



(4.50) 



4.4 Guessing the exact g(r) 



Let us try to gain more information about the power-law behavior in (4.1), that is about the 
expansion in 2— D of the correction-to- scaling exponent cu. Power-law behavior forces the 
series (4.50) to turn into some exponentially decaying function g(r) as can be seen from the 
asymptotic form of g(z): 



In order to check the latter equation note that 

poo 

flM z ) = zl+f / dr ex Ph 2 e ' r ~ (l+w/e)H = r(l+-) + 0(e~ z ) 
Jo e 

dr exp[-z e _r - (l+u/e)r] + 0(e~ z ) 
-u/e r) T 



(4.51) 



-u>/e 



r(i+ 



£ ' n=0 
1 X 



dr 



{-Uj/e) n 2 (n+l) 



r(i+ 



n! 



ii 



■ exp[— z e r — r] + 0(e" 



(4.52) 



£ ' n=0 

where it is understood that uo is expanded in powers of 2 — D. 



Let us now test a possible form of the exact g(r). It should satisfy the following properties: 

(i) In the limit of D = 2 the exact result r d / 2 /T(^) emerges. 

(ii) For D < 2 the corresponding g(z) has a finite fixed-point value together with a strong 
coupling expansion. Especially, the ansatz should interpolate to the limit D—l, which 
corresponds to a Gaussian polymer closed to form a ring. The strong coupling expansion 
of the renormalized coupling of a closed chain interacting with a 5— potential is easily 
obtained from the factorizability of loop integrals in D = 1 (see for instance [33]). The 
result is: 

°° / 1 \ n 1 

g{z)=e 1 + 



71=1 



T(e)zJ r(l-ne) 



(4.53) 
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(iii) It is consistent with the expansion (4.47). 
The (non-unique) ansatz is 



l-S(D,r) e-7' 



ujje 



d/2 



g(r) = C 

where S(D,r) is analytic in D = 2 of the form 

CO 

S(D,r) = l + "rJ2Sn(r)(2-D) 
and each S n (r) has a Laurent expansion 



(4.54) 



(4.55) 



n=l 



^max 



(4.56) 



3 — ^min 



Note, that in the limit of D — > 2, the expression (4.54) gives Cr d l 2 , while for D < 2 it yields 
upon integration the form (4.1), ensuring both properties (i) and (ii). Let us finally check con- 
sistency with the expansion (4.47) up to the second order in 2— D: Inserting 



u/e = u 2 (2- Df+O '(2- Df 



(4.57) 



(the linear term in (2—D) has to vanish 2 ) into the ansatz (4.54) and expanding to second order 
in 2—D provides 



~ g {r) = C r d/2 



l--(5i(r)(2-£>)4 



d-2 



u 2 



S 1 (r) 2 + S 2 (r) ){2-D) 2 + 



(4.58) 



Explicitly, (4.47) becomes to second order in 2—D 

g(z) = f* +2 {z) - (2 - D) d - \Cfl+\z) - Cf d (z)} 

+ (2 - Df d - [2Cf ft + \z) + (C 2 - 4Cf)/ 1 rf+2 (^) + (3Cf - C 2 ) f d (z) - Cf ft\z) 



A 2 



+ (2 - D) 2 j [2C 2 C f d+ \z) - (2C 2 C + C 2 )ff +2 {z) + 2C Z ff{z) - C' ff 2 (z) 
+ 0(2 - Df . (4.59) 

From this, the first coefficients of the (2— _D)-expansion of g(r) are obtained. They read 



^r2 



g{r) 



-d/2 



d 



1 + (2-D)-C 1-— - (2-D) 



2r 



^C!r+^(c 2 -4C!)-|(2C!+C 2 ) 



+ 



f(-^)4tf) r - 1 4(|- 1 )(e f+ ^)^]},4.60) 



2 This is due to the fact that the order (2 — D) term in g(z) scales identically in z as the leading term. Only the 
order (2 — D) 2 diverges more strongly. 
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Comparing (4.58) and (4.60), one identifies C = 1/F(^), Si = -C(l-f and u 2 = 2C 2 C , 

where C c (x):—C(x) — C. Note that the terms proportional to C 2 in S 2 (r) mostly cancel with 
Si(r) 2 , a sign that the ansatz catches some structure. 

The diagrams to be calculated at this order are C and C 2 . (see appendix B). On a manifold of 
toroidal shape, which is equivalent to periodic boundary conditions, two discrete sums have to 
be evaluated: 



c= Sd 



4vr 2 



. ^ 1 2tt 
^ p ~ (2 - D) 



= -0.44956 + 0.3583 (2 - D) + 0(2 - D) 2 (4.61) 



C c = j!h E ^ = 0.152661 + 0(2-/}). (4.62) 

fceZ^O 

With the results given above, this leads to 

to = 2eCf(2 - D) 2 + 0(2 - D) 3 = 0.305322 e (2 - D) 2 + 0(2 - Df , (4.63) 

which can be compared to the exact result for D = 1 (polymers): oo = e. As a caveat, note that 
the above scheme is not unambiguous in the sense that the second order term proportional to r 
in (4.61) could in principle either be attributed to cu 2 or S 2 . However, any ansatz in (4.54) will 
provide an uj, whose expansion starts at least quadratically in 2 — D. Though (4.54) is the best 
ansatz that could yet be found ensuring properties (i)-(iii), the precise form of constraints on the 
scaling function S remains to be discussed in order to settle this question. 



5 Conclusion 

In this work we refined the analysis of a D-dimensional elastic manifold interacting by some 
5-potential with a fixed point in embedding space. Starting from the perturbation expansion 
of the effective coupling of the problem, in a first step, we performed a new calculation using 
a modified regularization prescription: Evaluating loop integrals in fixed space dimension on 
a manifold of finite size enforced the introduction of a microscopic cutoff as soon as Z) = 2. 
This way, we recovered the complete summability of the perturbation theory in this limit and 
confirmed the strong coupling behavior as found previously in an analytic continuation from 
below D = 2. In the strong coupling limit, corresponding to strong repulsion or equivalently to 
large membrane sizes, the effective coupling diverges logarithmically as a function of the bare 
coupling z yielding a vanishing correction-to- scaling exponents. Analyzing the RG /3-function 
we found that it tends to zero at infinite bare coupling z as < d < 2. The renormalization 
group flow then tends to a fixed point, and the theory becomes scale invariant in this limit. Due 
to the logarithmic divergence of the effective coupling, however, the corresponding zero of the 
/3-function in terms of the latter is, too, shifted to infinity. This is a quite remarkable result 
showing that the scaling behavior of the system is accessible only to an all order treatment and 
deviates qualitatively from any finite loop expansion, be it within a minimal subtraction scheme 
or at finite e. Especially, the logarithmic growth of the effective coupling signals the limiting 
behavior of a scale-invariant theory. 

The result in D = 2 is completely independent of the regularization procedure. This does no 
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longer hold true beyond the leading order, which should be accessible to an expansion in 2 — D. 
We constructed its first order in a specific regularization scheme in [33]. While this reproduces 
qualitatively correctly the known result in D = 1, it suffers from a renormalization scheme, 
which neglects the boundaries of a finite manifold. We used a hard cutoff in position space, 
while working with the infinite D-space correlator. It seems that only in an e-expansion this 
procedure is systematic. 

Now, in a second step of we have overcome this problem by constructing the 2 — D-expansion 
on a manifold of toroidal shape of finite size, thus imposing periodic boundary conditions on 
the field. There is no further infrared cutoff necessary. We have carried out the expansion of 
the renormalized coupling up to fourth order in 2 — D, revealing the general structure of the 
expansion. It is important to point out that in considering g as a function of the bare coupling, 
the limits D — > 2 and strong coupling (z — > oo) can not be interchanged. While g tends to 
infinity as z does in D = 2, we expect finiteness of this limit as soon as 2 — D > and the 
existence of a strong coupling expansion as found for polymers (D = 1). We were able to 
guess an exact g(D, z) as a function of z and the internal dimension D, which satisfies these 
properties and which can be reconciled with the available expansion in 2 — D by an appropriate 
matching of its free parameters. Though it turned out that due to an ambiguity in the matching 
of parameters the precise power-law behavior of the effective coupling below D = 2 can not 
yet be isolated, we found that for closed manifolds the expansion of uj in powers of 2 — D starts 
at least quadratically as D < 2. 

The exponent is closely related to observables, which are accessible through Monte-Carlo ex- 
periments. These are for instance plaquettes-density functions at the repelling potential on a 
membrane avoiding a single point. 

While results for the pinning problem are interesting on its own, the main motivation is certainly 
to obtain a better understanding of self-avoiding polymerized membranes. Preliminary studies 
[40] indicate that this problem can also be attacked by the methods developed in this work. 
This would be very welcome to settle discrepancies between field theoretic results on one side 
[17,16,41] and numerical results (e.g. [32]) on the other. 



A The propagator 



The regularized difference correlator is defined as 



C a (x) = G a (0) - G a (x) , 



(A.1) 



where G a (x) denotes the usual two-point correlator, which is obtained from 3 : 




(A.2) 



3 Strictly speaking, we have to consider the propagator on the torus, as is done in appendix B. However this 
does not make any difference for the purpose of our argument. 
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Here, short-wavelength modes are suppressed through a soft cutoff procedure. Introducing a 
Schwinger parameterization for the evaluation of the integral in (A.2), 



l/a 2 

= f ds 8 D/2-l e -s4 

(2V^) D J 



(A.3) 



where s = l/(t + a 2 ), we obtain for (A.l): 

C a (,) = ^/ ds ^(l-e^). 



(A.4) 



Further evaluation leads to: 
(i) D = 2: 



1 / 2 2 \ -i 



C . W T + r(0,^) + ln- ^-ln-. (A.5) 



(ii) D < 2: 



WW (2 _^ 2 7rD/2 ^ (2 _ D) 2 d - 1 tt^/ 2 ° (2-D) 2 D ~ 1 n D / 2 (2-D) 2 n D / 2 

D^2 



ir^oo \X\ tt D^2 J- , X 



S D (2-D) (2-D)2 d - 1 tt d / 2 2tt a 

T(z, a) denotes the incomplete T-function: 



In - . (A.6) 



F(z,a) = I dtt* _1 e _ * . (A.7) 

Especially: 



as long as D < 2. 



B Calculation of the diagrams in the (2 - L>)-expansion 

In this section we calculate the diagrams which appear in the 2 — D expansion on the torus of 
size L = 1. It turns out that to obtain C and C 2 we need to evaluate two sums over discrete 
wave-vectors due to periodic boundary conditions on the torus. Let us first derive the latter 
before turning to the explicit evaluation. Starting from the definition of the difference correlator 

C(x), 

C(x) := G(x) - G(0) , (B.l) 
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where G(x) is the usual two-point correlator, we obtain C(x) through an inverse discrete 
Fourier-transformation from G(k) = 1/k 2 , which reads: 



C(x) = \ ( x - el& ) > k = 2nn, neZx Z\{0} . 

k^O ^ 



(B.2) 



Performing the averaging procedure 



C(x) = [ C{x) , 

J x 



(B.3) 



where J x 6 lkx = 8? is to be taken into account, the calculation of C(x) reduces to 

CR = /i:=Vi, k = 2nn, 
r - ^ k 2 



(B.4) 



where k is D-dimensional, and the indices are integer and running from — oo to oo, n = 
being excluded from the summation. Of course, in the expansion in powers of 2 — D we need 
an analytic continuation to real values of D. Finally, to obtain C(x) we have to subtract C^(x) 
from C(x). Due to our normalizations: 



C(x) = S D C(x) 



2tt(2 -£>)/' 



(B.5) 



where So denotes the volume of the unit sphere and C(°)(x) = 1. 
Turning to C 2 (x), we first note that within our normalizations we have 



S» 2 C 2 (x) = (C(x) - C(°)(x)/(27r(2- J D))) 2 = C(x) 2 - 2 



C(x) 



2w(2-D) (2ir(2-D))< 



(B.6) 



and 



Sp 2 C(x) =C(x) -2- 



C(x) 



+ 



2n(2-D) (2vr(2- J D)) 2 
according to (B.5), such that 

Cp) = C 2 ^) - £{xf = S 2 D (cHx) - C{x) 2 ) ■ 
Knowing already the sum to be evaluated to obtain C, (B.4), what is left is: 
CHx) = Jc 2 {x) = Jj2 J: if" - l) (e«* - l) 



(B.7) 



(B.8) 



k^O P^O fc^o 



^ k 2 

k^O 



(B.9) 
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Therefore, 



S D 2 Cl(x) = I 2 :=J2i- 4 > h = 27rn t . 

-» rv 



(B.10) 



Let us first calculate Ii. Introducing a Schwinger parameterization we have: 



1 = (27r) 2 ^2 = E 



7t j — — OO 



n 2 (2ttY 



rij — — oo 



ds e~ 



(2tt) s 



ds 



E 



(B.ll) 



where it is to be noted that the sum in the last line is only one-dimensional. Furthermore, from 

now on it is clear, how I\ is analytically continued to real values of D. 

In order to evaluate this sum, we will make use of a Poisson-transformation, which reads: 



E 



-A(n-z/2f 



_ \ — +fwlz 

A ^ 



(B.12) 



l=— oo 



The contribution from / = is the approximation of the l.h.s. through a Gaussian integral. 
Our aim is to calculate the coefficients of the 2 — D expansion of Ii numerically using some 
algebraic manipulation program. Then, the integration interval in (B.ll) has to be made finite. 
This is done as follows: For any s > we have 



h 



»0 



(27T) 



ds 



ds 



D 



E 



E 



D 



+ 



+ 



(2 



oo 



ds 



D 



E 



so 



(27T) 



ds 



D 



e- 2/s -1 



(B.13) 



For any finite s > 0, the sum in the r.h.s. integral can be truncated at some finite n max for 
all s G [0, s ]. For the first integral (corresponding to small values of s) we make use of the 
poissonian formula (B.12) with z = 0: 



E - , E ^" 



2 /s 



(B.14) 



l=— oo 



Inserting this into (B.13), the sum in the first integral can be truncated at some finite I as well, 
such that one may approximately write: 



h 



ds 



D 



e-^ 5 -1 



l = -ln 



SO 



+ 



(27T) 



ds 

^2 



D 



e- 2 /M-l 



vn=— n m ax 



(B.15) 
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Choosing s in a way that l max can be set equal to zero the l.h.s. integral can be evaluated 
analytically: 



h 



2tt d / 2 
2-D S ° 



D/2-1 



- Sn 



SO 



+ 



(27T) 



ds 



D 



e~" 2 /M-l 



(B.16) 



There is a pole in 2 — D, which can be easily subtracted expanding the expression in powers of 
2 — D. The pole is 

/ ' = M2^) +0 « 2 - D »°»- (B ' 17) 

The precision of the machine that we used to evaluate (B.16) was sufficient in a way that we 
could select s from an interval, such that the sum appearing in the integrand could be truncated 
at some finite n max and the result was independent from the precise value of s within the desired 
order of accuracy, therefore, justifying the approximation in (B.15). Setting for instance s = 
1.9 and n m „ = 20 we obtain with Mathematica®: 



h = 



2tt(2 - D) 



0.715497(1) - 0.00457046(1)(2 — D) + 0((2 - D) 



(B.18) 



On the torus we scaled the square root of the volume of the D-dimensional unitsphere into the 
field. Accordingly, comparing with (B.4) and (B.5) we then find: 



C = -0.44956(1) + 0.3583(1)(2 - D) + 0((2 - D) 



(B.19) 



Let us turn to the evaluation of I 2 following the same strategy as above. Again, setting L = 1 
and introducing a Schwinger parameterization leads to: 



h = 



(27T) 4 
1 



ds s 



ds s 



D 



E 



D 



7 E 



-7T 2 1 2 /S 



l=-ln 



+ 



ds 



D 



E 



-n 2 /s 



(B.20) 



where we have once again applied the Poisson-transformation (B.12) with z = on one part of 
the integration interval and truncated both series at some finite values n max and Z max . 
There is no pole in 2 — D. Since I 2 appears at second order in 2 — D we only need its value at 
D = 2. s has to be chosen from an appropriate interval. Setting n max = n max = 10 and s = 1.1 
we obtain with Mathematica®: 



h = 0.00386695(1) + 0((2 - D)) , 



(B.21) 



or, due to the rescaling by S 2 D , 



C2 = 0.152661(1) + 0((2-D)) 



(B.22) 
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